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Abstract 



Spin fluctuations enter the calculation of the superconducting transition 
temperature Tc only in the next-to-leading order (i.e., in 0{1/N'^) of the \/N 
expansion of the t — J model. We have calculated these terms and show that 
they have only little influence on the value for Tc obtained in the leading 
order 0{1/N) in the optimal and overdoped region, i.e., for dopings larger 
than the instability towards a flux phase. This result disagrees with recent 
spin-fluctuation mediated pairing theories. The discrepancies can be traced 
back to the fact that in our case the coupling between electrons and spins is 
determined by the t-J model and not adjusted and that the spin susceptibility 
is rather broad and structureless and not strongly peaked at low energies as 
in spin-fluctuation models. 

Relating Tc and transport we show that the effective interactions in the 
particle -particle and particle-hole channels are not simply related within the 
1 /N expansion by different Fermi surface averages of the same interaction as 
in the case of phonons or spin fluctuations. As a result, we find that large 
values for Tc and rather small scattering rates in the normal state as found 
in the experiments can easily be reconciled with each other. We also show 
that correlation effects heavily suppress transport relaxation rates relative to 
quasiparticle relaxation rates in the case of phonons but not in the case of 
spin fluctuations. 

PACS numbers: 74.20-z, 74.20.Mn, 74.25.Ha 



I. INTRODUCTION 



Calculations of the superconducting transition temperature Tc and other specific prop- 
erties of the cuprates are usually based on the Hubbard model. One approach assumes 
that the effective Hubbard repulsion U is smaller or comparable to the electronic hopping 
ampliutde t so that approximations which become exact at weak coupling are taken as a 
basis for calculations of properties of the cuprates. The FLEX approximationB, which sums 
bubble and ladder diagrams in a self-consistent way, is such an approximation which has 
been rather successful in reproducing properties of the cupratesi. A rather different and 
certainly more realistic view of the cuprates assumes that U is substanially larger than t 
at all relevant dopings. As a result the strong-coupling limit of the Hubbard model applies 
which is equivalent, except for three-center potentials, to the t — J model. Cuprates are 
then no longer treated as moderately correlated metals as in the FLEX approximation but 
as doped Mott insulators. Formally this means that self-energies, effective interactions, etc., 
are not expanded in powers of U but in 1/f/ and that, as a new point, the constraint of 
having no double occupancies of sites has to be taken into account. Throughout our paper 
we adopt this second point of view and thus assume that the strong-coupling limit applies 
to the cuprates over the whole doping region including also the overdoped case. 

The t — J model has no small parameter which would allow a systematic perturbation 
expansion. One way out is to extend this model in such a way that a small parameter arises 
and that the properties of the original model are obtained as power expansion in this small 
paramter. One convenient generalization consists of replacing the two spin components of 
the t — J model by N components and, correspondingly, its orginal SU{2) symmetry by the 
symplectic symmetry group Sp{N/2). In this way properties of the physical N = 2 model can 
be expanded around the N = oo limit of the extended model which describes renormalized, 
but non-interacting, particles. The interaction between the particles can then be taken 
into account by means of a systematic expansion. We have formulated this expansion 
directly in terms of X-operators which respect the constraint via their commutator and 
anticommutator rules. Applying this method to superconducting instabilities we foundlHi 
for a square lattice relevant instabilities only in the d-wave channel {T^ representation of 
the point group C^^ of the square lattice) due to the Heisenberg term. Putting N = 2, 
the corresponding transition temperature Tc was about 0.03 in units of the nearest-neighbor 
hopping element t for a doping value away from half-filling of 5 = 0.17. 

The above calculation took into account only the leading order of the expansion. 
Thus the anomalous self-energy was calculated in 0{1/N), or equivalently, the kernel B of 
the linearized gap equation in Oil), since an overall factor 1/A^ has been taken out in the 
gap equation (see Eq.(40) of Ref.Q, or Eq.(l) below). Although this order contains also spin- 
fiip contributions it does not contain the well-known RPA terms for the spin susceptibility. 
These terms are, of course, present in the 1/A^ expansion but their leading contribution (one 
bubble) yields only 1/A^ terms to and thus has been dropped in Ref.i. 

Presently many experimental data are interpreted as giving evidence for a spin- 
fluctuation mediated pairing mechanism for high-T^ superconductivityH. It is thus in- 
teresting to include the above terms, at least, if they are of 0{1/N) in 0, in the calculation 
of Tc. It is the aim of the present communication to do this and to consider their effects 
on Tc and the resistivity. Using the 1/N expansion of the t — J model these terms can 
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easily be identified. The corresponding coupling strength between electrons and spins is 
determined uniquely by the model parameters and has not to be chosen in an ad hoc way. 
At the same time such a calculation yields a check for the validity of the expansion 
for T^: If, for instance, the next-to-leading terms are of the same magnitude as the leading 
ones, after putting = 2, a rapid convergence of the series cannot be expected. Of course, 
the comparison of the two orders would be most convincing if all the 0{1/N) terms of 
could be taken into account. The resulting expressions are, however, rather involved and 
their evaluation has to be postponed. In the following we will restrict ourselves to that part 
of B, called 0^, which is of 0{1/N) and describes RPA-like spin fluctuations. These terms 
are derived in section 2. Section 3 discusses properties of the resulting spin susceptibilty and 
compares them with a different theoretical approach and also with experimental data. In 
section 4 properties of G<j, which are relevant for superconductivity, are discussed and their 
influence on Tc is calculated. In section 5 the consequences of 6 for transport properties 
are investigated. Section 6 contains our conclusions. Throughout the paper we will restrict 
ourselves to dopings 5 > 5^, where 5c denotes the critical value for the instability towards a 
flux or bond-order phaseSli0'0. 5c is about 0.15 for J/t = 0.3. It has been arguedS that 5c 
corresponds to optimal doping so that we deal in the following with optimal doping and the 
overdoped regime. 



II. BASIC EQUATIONS 

Using a 1/A^ expansion for the t— J model the linearized equation for the superconducting 
gap Y^an can be written aM 

s»W = -^i:e(*,*');3j:^E.„(.'). (1) 

A^ denotes the number of electronic degrees of freedom per site. We assume that they consist 
of A^/2 replica of one orbital with two spin components and A^ can be assumed to be an 
even integer for convenience. Nc is the number of primitive cells and k the supervector 
k = (co'„,k), where uJn denotes the fermionic Matsubara frequency uJn = {'2n + l)vrT, k the 
wave vector, and T the temperature. e(k) is the one-particle energy 

e(k) = ^t(k) - J(k) • i- X: cos{p.)f{e{p - /i)), (2) 

with 

i(k) = -2\t\{cos{k^) + cos{ky)), (3) 

J(k) = 2J{cos{k^) + cos{ky)). (4) 

/X is the chemical potential, / the Fermi function, 5 the doping away from half filling and the 
lattice constant of the square lattice has been put to one. Note also that the usual definition 
of the hopping matrix element in A^ = 2 theories corresponds to \t\/2 in our notation. 
Generally accepted values for \t\/2 and J are about OAeV and 0.15eV, respectively. In the 
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following we will use \t\/2 as energy unit. The kernel represents the effective interaction 
in the particle-particle channel. Its leading 0(1) contributions consist of four different terms 
0(1)^ 0(2)^ 0(3)^ 0(4)^ gi^gj^ explicitly in Eqs.(42)-(52) of Ref.l. The first one describes 
the instantaneous part, the second one fluctuations in the charge density, and the third and 
fourth ones mixed spin and charge fluctuations due to the anomalous part of the vertex. 

0(2) 

originates from single particle and collective electron-hole exitations in the charge 
channel. A similar contribution 6^ is also obtained in the spin channel 

Qs{K k') = -(t(k') + J(k - k'))7s(fc, k' -k), (5) 

where the sign for singlet pairing has been chosen. Expanding 7^ in powers of the 
leading 0(1) contribution is just a product of delta-functions. This contribution has already 
been taken into account in The 0(1/A^) contribution of 7^ contains the term 

7^(^:1,^2) = ^(t(kl) + J(k2))x(°nfc2), (6) 

with 

^ e(k + q) - e(k) - lujn 

Inserting Eq.(|D into Eq.(^ we have the following additional 0(1/A^) contribution to the 
kernel 9^ due to spin fluctuations 

e.(fci, k2) = -^(t(k2) + J(ki - k2))(t(ki) + J(ki - k2))x^°nfci - k2). (8) 

Both the susceptibility x^^^ cind the coupling functions on the right-hand side of Eq. (|D 
are different from what one may expect from a simple spin-fluctuation mediated theory for 
superconductivity. The coupling functions (first and second factors in the above expression) 
contain not only the potential J but also the kinetic parameter t. The origin of the latter 
lies in the fact that fermionic Hubbard operators are not usual creation and annihilation 
operators and generate many-body interactions in the perturbation expansion. Secondly, 
one would expect the RPA expression for the spin susceptibility on the right-hand side of 
Eq.(|D and not the equivalent of just one bubble, as described by Eq.(^. The reason for 
this is that in RPA-type diagrams each bubble is multiplied by J/N , i.e., entering as an 
0{1/N) contribution. Thus the infinite summation over bubbles in the RPA is not in the 
spirit of the 1 /N expansion: Before the product of even two bubbles can be included one has 
to sum up infinite series of diagrams of different types. This can be seen from the 0(1/N) 
expression for the spin vertex 7^ in Refs.lillii. Since already exhausts more than half of 
the exact sum rule for the spin susceptibility (see below) we retain from the 0(1/N) terms 
to only the term 6^, as given by Eq.(RD. 



III. SUM RULE AND PROPERTIES OF THE SPIN SUSCEPTIBILITY 

For a general N the spin operator at site i is replaced by 
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Sr - S^u), (9) 

where z/ is a flavor index counting the N/2 copies of the original orbital and a a Cartesian 
index counting the three vector components of the spin. Expressed by Hubbard operators 
the z-component of the spin operator acting on the copy u reads 

^f(z.) = i(xr'^'^-xr'^^). (10) 



X-^k,u) = -{-z) Y: I dte^-'{[S^{t,u),S^_^{0,u)]), (11) 



The spin susceptibility x""^ is defined in Fourier space for a general N by 

u=l...N/2' 

where (..) denotes a thermodynamic average and [...] the commutator. Inserting Eg . (p!0D into 
Eg. (p!TD , using the fluctuation-dissipation theorem, summing over k and integrating over lj 
we obtain for the z-component of the imaginary part of x iii the normal state 
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-^E/ Imx-ik,co)il + n{co)) = ^ ^ ((Xf'^'^ - X^'^T)- (12) 

k -^-"^ u=l...N/2 

n{uj) is the Bose function. The right-hand side of Eg.(|l^) can be evaluated exactly for the 
two cases N = 2 and N = oo. For N = 2, which corresponds to the physical case of one 
orbital with two spin components, the X-operators are projection operators, so that Eg.(|l2D 
becomes, using also the constraint X^^ + X^^ + XP = 1, 

-TtE/ dujImr\\,,u){l + n{uj)) = -{l-5), (13) 

where 5 is the doping away from half-fllling. Eg.(0) is exact and generalizes the spin sum 
rule to a general doping 5. Experimental data are usually given in terms of the magnetic 
susceptibiltiy which is the product of the spin susceptibility and fi'^/i^, where g and iib 
are the gyromagnetic ratio and the Bohr magneton, respectively. Following a widely used 
notation we denote the negative imaginary part of the magnetic susceptibility, analytically 
continued to the real axis iun ^ u + ir], by x"(k, cu). Inserting g = 2 the sum rule Eg . (pTsD 
assumes then the form 

— / rfu;x"(k,u;)(l + n(^))=/i|7r(l-5). (14) 

For N > 2 the projection properties of the X operators are lost and non-diagonal averages 
on the right-hand side of Eg.(|l^) are nonzero. However, at large X's the non-diagonal 
averages are by a factor 1/X smaller than the diagonal ones and thus can be omitted. The 
diagonal ones can be calculated using Eg.(|^) and one obtains intead of Eg.(p^) 
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— Y du;x"0^,u;){l + n{u;))=^,ln{l-6')/2. (15) 

At large dopings Egs.(|l^) and (|15D agree approximately but at zero doping the right-hand 
side of Eg . ([15|) is only 1/2 of that of Eg . ([l^) . The right-hand side of the above sum rules, 
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considered as a function of A^, interpolates between the above two extreme cases N = 2 and 
N = oo. This means that terms ~ 1/iV and of higher orders contribute to the sum rule and, 
if evaluated at = 2, cancel the factor (1 + 5)/2 on the right-hand side of Eq.(]l5D. This 
implies, in particular, that the sum rule for the spin susceptibility depends on the strength 
of the interaction between the particles. 

Fig. 1 shows theoretical and experimental data for the momentum averaged magnetic 
susceptibility x"{^)y defined by 

X» = ^Ex"(k,c.). (16) 

^ k 

The dot-dashed line represents the susceptibility at large A^'s, x'Li^)- The dashed curve in 
Fig. 1 is a plot of the function Xp{^) which is associated with the following phenomenological 
expression for the spin susceptibility in high-Tc compounds!^. 

Here we have used the values listed in Table II of Ref.li for the various parameters for 
6 = 0.2 (second-last column in Table II). The open circles in Fig. 1 are experimental 
points for the momentum averaged susceptibility divided by 2 in underdoped YBaCu^OQ,^ 
at T = 80K. Analogous experimental data have been presented in Refs.llHil for various 
dopings in YBa2Cu30j^x- The data for YBa2Cu30Qj (Fig. 6 of Ref.0), covering the 
energy window between and 70 meV, are smaller by roughly a factor 2, but otherwise 
rather similar: They also show a first maximum at around 30 meV and a second maximum, 
which lies at a somewhat lower energy of about 60 meV and is broader and less pronounced 
compared to the experimental data shown in Fig. 1. Increasing the doping in YBa2Cu307-x 
spin scattering in the normal state becomes weaker and finally undetectable in the optimally 
doped case0. x'L exhausts a little more than half of the sum rule Eg. ([T^) . It spreads in 
frequency over the bandwidth of about 0.7 eV and exhibits a broad maximum near 0.3 
eV. Integrating x'Li^ + n{u!)) up to the largest measured frequency the dot-dashed and the 
experimental curve exhaust about 1/6 of the sum ru However, the experimental points 
for the oxygen contents 6.6 and 6.7 show a dramatic shift of spectral weight from high to low 
frequencies which is absent in x'^- The dashed curve in Fig. 1 is of similar magnitude as the 
experimental curve between and about 70 meV. However, it misses the strong decrease of 
the experimental curve above 70 meV and is therefore much too large at high frequencies. 

Fig. 2 shows data for ^"(Q,^^) at Q = (tt, vr) where the susceptibilities assume either 
their maximal values in k- space or are near their maximal values in case where a splitting 
into incommensurate peaks positioned near Q occurs. The filled circles are experimental 
values from Ref.0 for Y Ba2CusOQ,7 at T = 200K. The dot-dashed and dashed curves in 
Fig. 2 represent x'LiQy^)^ multiplied by 10, and Xp{Q,^), respectively. Note that there is 
a scale change of about a factor 50 along the y-axis between Figs. 1 and 2. This means that 
the experimental spectral weight of the spin susceptibility is strongly concentrated near Q 
and at low energies around 20 or 30 meV. This strong localization in momentum and energy 
is absent to a large extent in x'L and to a minor degree also in Xp- Plotted as a function of 
momentum for a fixed frequency x'L shows well-pronounced incommensurate peaks near Q 
in spite of the fact that x'4)(Q) ^) as a function of uj is according to Fig. 2 rather small and 
structureless. 



6 



The showing of experimental points for the underdoped regime in Figs. 1 and 2 does 
not mean that ^'4 should reproduce these experimental spectra. As mentionned in the 
introduction our expansion holds in this form only for dopings 6 larger than 6c where the 
instability towards a flux phase occurs, i.e., as discussed in Ref.0, in the optimally doped and 
overdoped region. Presently spin scattering in optimally doped Y Ba2Cu^0j is too weak to 
be detectable by inelastic neutron scattering with the instrumentation currently availablelll. 
Recently, an upper limit for x"(Q,cij) of Abfi^/eV has been given for the instrumentally 
accessible energy range 20meV < u < 50meV^ which is substantially lower than the 
dashed line in Fig. 2. This suggests that a low-frequency peak in x" is present only at low 
dopings and that x" optimally doped and overdoped samples extends over a large energy 
region because of the exact sum rule Eg. ([131) . One thus may argue that x'L represent 
well the spin susceptibility in the optimal and overdoped regimes and that the experimental 
points as well as Xp Figs.l and 2 only apply to the underdoped case. 

Another quantity characterizing x is the coherence length C,- In an isotropic two- 
dimensional description the decay of the the real part of xi^^^ = 0) i^i real space is given 
by Ko{\r\/^), where Kq is a modified Bessel function of zeroth order. The dashed and dot- 
dashed lines in Fig. 3 represent this Bessel function for ^ = 2.3 and ^ = 1.0, respectively. 
The diamonds are the absolute values for the real part using xp at various neighbor sites 
of the square lattice. The diamonds lie perfectly on the Bessel function with ^ = 2.3 as 
expected. The circles in Fig. 3 have been calculated with the large-N susceptibility Xoo- 
Because of the anisotropy of the lattice not all of the calculated points lie on the curve for 
the Bessel function. Fig. 3 nevertheless shows that the corresponding ^ is near the value 1. 
Repeating the same calculation at finite temperatures we find that ^ varies with tempera- 
ture on the scale of t and thus is independent of temperature in the range of experimental 
interesting temperatures. Experimental values for ^ from inelastic neutron scattering are 
independent of temperature and have been summarized in Fig. 2 of Ref. clS db function of 
doping. For optimally doped YBa2Cu30j-x C is only little larger than 1 which also agrees 
with NMR dataS'i. This supports the idea that Xoo represents a sensible approximation 
for the spin susceptibility in the optimal and overdoped regime. 



IV. PROPERTIES OF Os NEAR THE FERMI SURFACE AND ITS INFLUENCE 
ON THE TRANSITION TEMPERATURE Tc 

A relevant quantitiy for superconductivity is G(k, k',cij = 0) with momenta k and k' 
put to the Fermi surface. Fig. 4 shows this quantitiy for a fixed first momentum at k = 
(2.465, 0.309) as a function of the second momentum k'. k' moves counterclockwise around 
the Fermi line passing through the points X,Y, {X,Y) along the positive (negative) x- and 
y-axis, respectively. The dashed line corresponds to the kernel g'^Xp/^: using the coupling 
constant g = O.QeV from Table II of Ref.i. The dot-dashed and solid lines in Fig. 4 describe 
the kernel and the full 0(1) kernel B calculated in Ref.i, respectively. The doping is 
always 6 = 0.17 and J equal to 0.3. The additional factor 1/A^ in B^ has been taken into 
account using the physical value N = 2. The curves are not symmetric with respect of X 
or X because the first argument k has been put to a general point somewhat away from X. 

First we note that all three curves have maxima at Y and Y, i.e., for momentum transfers 
of Q = (tt, tt) or of — Q. The dashed and dot-dashed curves correspond to purely repulsive. 
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the solid one to predominantly repulsive interactions. As a result d-wave pairing is preferred 
in all three cases. The variation in the solid curve along the Fermi surface, which for 
d-wave pairing is the relevant quantitiy, is about two to three times larger than for the dot- 
dashed curve. This means that the spin-fluctuation mediated interaction is by this factor 
smaller compared to the sum of all 0(1) contributions to the effective interaction. This 
already indicates that the 1/N expansion makes sense insofar that the terms decrease with 
increasing order in The dot-dashed line in Fig. 4 exhibits some fine structure near 

the point Y which is caused by the splitting of the peak at Q into incommensurate peaks 
displaced slightly away from Q. The appearance of these incommensurate peaks due to the 
Fermi-surface geometry has been discussed previously in Ref.^. 

The kernel 0(k, k', 0) represents the effective interaction in the particle-particle channel 
at zero frequency. The solid line in Fig. 4 thus can be compared with the upper curve in 
Fig. 8 of Ref.S which has been obtained from Monte Carlo simulations. Though both kind 
of curves show minima at the X and maximma at the Y points the energy scale of the curve 
m Ref.@ is one order of magnitude larger than in our case and also has a large momentum 
independentoart. These differences can roughly be explained by the fact that the energy 
scale in Ref.c3 is set by the Hubbard constant U and in our case by J(k). The parameter 
values used in Ref.@ thus correspond more to the weak-coupling case whereas we deal 
with strong-coupling. Fig. 4 also shows that the momentum dependence of the effective 
interaction is similar to that of the spin susceptibility which determines the momentum 
dependence of the dashed and dot-dashed curves. However, one cannot conclude from this 
that the effective interaction is caused by antiferromagnetic spin fluctuations because the 
solid line does not include them. 

The dashed curve in Fig. 4 is at least one order of magnitude larger than the dot- 
dashed curve, especially, near the point Y. Though both terms describe interactions due to 
spin fluctuations their strengths differ roughly by one order of magnitude. To understand 
why these curves are so different in magnitude we characterize the coupling in the large-N 
expression Eq.(|) by an effective coupling constant g. Since both, t and J, depend on the 
momentum we define g by the following average value appropriate for d-wave symmetry 

f =« e,(k,k',a; = 0) »3 / « x^°Hk,k',^ = 0) »3 . (18) 

<< ... >>3 denotes in the above equation a momentum average over the Fermi line with 
a weight function with d-wave symmetry. An appropriate weight function is given by the 
eigenvector associated with the lowest eigenvalue with d-wave symmetry. The explicit form of 
the corresponding kernel yielding the desired eigenvalues and eigenvectors has been described 
in Ref.i. This definition can be extended to any momentum-dependent operator. The Fermi 
surface average is then obtained as the lowest eigenvalue of the associated operator. The 
dot-dashed and the dotted lines in Fig. 5 show the numerator and the denominator of the 
right-hand side of Eq. ([I8| ) as a function of doping. One concludes from these curves that g 
depends only weakly on the doping and that it is about 0.5 to 1 eV in magnitude. If one 
puts t to zero in the numerator on the right-hand side of Eq.(|l^) the ratio is about 0.1 {eVy. 
This means that g is determined essentially by t and given approximately by the full band 
width. This result is a typical strong-coupling result because the effective coupling constant 
g is determined by the kinetic energy and not by potential interactions. One also concludes 
from this that the coupling constant g = 0.6eV used in Ref.S for the coupling between holes 
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and spin excitations is not too large and thus does not cause the difference in magnitude 
between the dashed and the dot-dashed curves in Fig. 4. Also shown in Fig. 5 are the 
lowest eigenvalues of d-wave symmetry for the total kernel of the 0(1) theory (solid line), 
the kernel (dot-dashed line, which is identical with the line representing the numerator 
on the right-hand side of Eq.(|18|)) and the kernel g'^XP (dashed line). The eigenvalues of 
the two spin-fluctuation mediated kernels differ by roughly one order of magnitude which 
reflects the ratio of the corresponding curves in Fig. 4. Using the more realistic valuS 
^ ~ 1 in xp the dashed curve would dramatically move upwards and be near the dotted or 
dot-dashed curves. Interesting is also that the eigenvalues of the leading order of the large-N 
theory (solid line) is only about a factor 2 larger than those of B^. 

Fig. 6 shows curves for the negative imaginary part of the kernel B as a function of 
frequency. This quantitiy is the analogue of the familiar function a'^F{uj) of Eliashberg 
theory for conventional superconductivity. The dependences on the momenta k and k' have 
been averaged over the Fermi surface using the lowest eigenvalue of d-wave symmetry of 
the appropriate operator. The solid line is the result of the large-N theory, omitting the 
instantaneous, attractive term. The dot-dashed and the dashed lines correspond to the 
spin-fluctuation mediated interactions described by Bg and g'^XP^ respectively. 

The d-wave average << ... >>3 over the Fermi line and the momentum average over 
the Brillouin zone may yield quite different results as can be seen by comparing the dashed 
and dot-dashed curves in Fig. 1 with the corresponding curves in Fig. 6. Density of states 
effect are vital in the average << ... >>3, giving momentum transfer q ~ Q a large weight. 
As a result the dashed and dot-dashed curves in Fig. 6 resemble more the corresponding 
curves in Fig. 2 than in Fig. 1. The dashed curve in Fig. 6 exhibits a strong maximum 
around 15 meV and a slow decay towards high energies with a typical decay constant of 
about 100 meV. Compared to the dashed curve the dot-dashed curve in Fig. 6 assumes 
everywhere only small values, is rather structureless and has a broad maximum around 0.5 
eV. Since the coupling constants g and g are of similar magnitude in the two cases the 
difference between the dashed and dot-dashed curves in Fig. 6 is due to the difference in 
Xp and Xooi especially for momenta q ~ Q. If one restricts oneself to energies between zero 
and 0.8 eV the integrated spectral weight of the solid curve is much larger than that of the 
dot-dashed curve. This implies that the next-to-leading contribution in the 1/A^ expansion 
is substantially smaller than the leading one which gives evidence for the usefulness of the 
1/A^ expansion. In the same energy interval the average spectral weight of the solid curve 
is substantially larger than that of the dashed curve. However, the solid curve becomes 
strongly negative between 1 and 2 eV. This negative part compensates to a large extent the 
positive spectral weight between and 0.8 eV if one calculates the real part of the effective 
interaction. As a result the retarded 0(1) contribution to B becomes rather small at low 
frequencies and the attractive instantaneous part causes mainly the superconductivity, as 
discussed in Ref.l. Such a compensation is absent in the two spin fluctuation curves because 
they are always positive and also very small at high energies. 

Using the numerical procedure described in Ref.i to solve the linearized gap equation 
Eq.(|l]) we find for Tc the values shown in Fig. 7. The circles describes Tc as a function 
of 5 for the 0(1) expression for B as discussed in Ref.i. If only the new term B<j is taken 
into account the corresponding Tc is given by the crosses. These Tc's are at least one order 
of magnitude smaller than the values for using the 0(1) expression. In accordance with 
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that the Tc's for all terms included in the gap equation (diamonds) are only slightly larger 
than those calculated solely from the 0(1) contribution. 

V. IMAGINARY PART OF THE SELF-ENERGY AND RELAXATION RATES 

The particle-particle and the particle-hole effective interactions are, except for a possible 
minus sign in the case of spin fluctuations, identical if the effective interaction is mediated by 
bosons with energies much smaller than the Fermi energy. Tc, quasi-particle life times and 
transport relaxation rates are then obtained by taking different Fermi surface averages of the 
effective interaction. The particle-particle and the particle-hole interactions in 0(1/N) of 
the t — J model are, however, different. Quasi-particle life times and scattering rates are not 
averages of the function G of Eq.(|I|) but of the corresponding function B^'' in the particle- 
hole channel. Following Refs.Ull the 0(1/N) contribution to the normal-state self-energy E' 
of the electronic Green's function G, defined in Eq.(lO) of Ref.i, is given by 

T.'{k) = ^Y.^^\k,k')g{k + k'), (19) 

with 

Qph ^ Qph ^ Qph^ ^20) 

ef (fc, k') = -7e(fc, k'){t{k + k') + J(-k')) - (t(k + k') + J(-k'))(t(k) + J{k'))xu{k') 

+ 7c{k,k') E Er{k)x2r{k')-lc{k,k') E X2r{k') E Xrs{k')E,{k) 
r=l r=l s=l 

- E Srik'k)xr.ik')xstik')Etik), (21) 

r,s,t=l 

with 

Sr{k', k) = -j2{l + x{k'))-^X2s{k'hc{k, k'), (22) 



Qtik.k') = go\k)Q'\t{k + k') + Jhk')){Lrj ) + NL{ )). (23) 

L is the Fourier transform of the equilibrium Green's function of two bosonic operators X{1), 
X(2), except for a minus sign, 

L(l, 2) = (T(X(1) - (X(1)))(X(2) - (X(2)))). (24) 

go{k) is the free quasi-particle Green's function {iun — e(k))~^ and Q the spectral weight 
(X°°(i))-|-(X^^(i)). At large N the second term in Q could be dropped because it contributes 
only higher orders in 1/N than those considered here. For the transition from large N to 
the physical value iV = 2 it is advisable, however, to keep also higher orders in 1/N in the 
inverse spectral weight in order to improve the convergence of the 1/N expansion at small 
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dopine. The remaining symbols in Eqs.(pTD, (p2D, and (^) have been defined in Eqs. (41-48) 
of Refl. 

The self-energy S' is associated with the physical Green's function G. Taking only 
in Eq. (pOD into account Eq. (|19D defines a different self-energy S related to the quasi-particle 
Green's function g defined in Eq.(15) of Ref.i. If k is near the Fermi surface (i.e., e(k) near 
the Fermi energy and the frequency iuJn small) Gf^ of Eg . (^31) is small due to the factor 
gQ^{k). This means that the two kind of self-energies are essentially equal near the Fermi 
surface. However they differ considerably far away fom the Fermi energy: The imaginary 
part of E' is there always negative definite as it should be in contrast to that of E which 
assumes there in general both signs. 

Fig. 8 shows various contributions to the negative imaginary part of E' as a function 
of energy for t = OAeV, J/t = 0.3, 6 = 0.17, calculated at very low temperatures. The 
momentum in the self-energy has always been averaged over the Fermi surface. The dotted 
line describes the pure spin-fiip contribution, i.e., the second term in Eq.(pID and the second 
L-contribution in Eq . (|23|) . The dot-dashed line in Fig. 8 contains the charge contribution 
to the self-energy, i.e., the first term in Eq.(pTD and the first L-contribution in Eq. (p^) . 
The solid line in Fig. 8 contains the contribution of all terms in 0^^, i.e., the full 0(1) 
contribution of E'. First we note that the sum of the dotted and dot-dashed lines in Fig. 8 
is not much different from the solid line. This means that the terms 3-5 in Eq. (|2T|) describing 
mixed charge and spin fluctuations yield only a small contribution to E'. The imaginary part 
of E' shows a double peak both for negative and positive frequencies. The peaks at lower 
frequencies are due to spin, the peaks at higher energies due to collective charge fluctuations. 
This is in agreement with the properties of the spin and charge excitation spectra of the t — J 
modelii. For relaxation rates and transport only the low-frequency part of the imaginary 
part of E' is relevant. Fig. 8 suggests that this part is mainly determined by spin and only 
to a lesser degree by charge fluctuations. 

Theories for high-Tc superconductors which are based on a retarded, boson-mediated 
pairing mechanism usually encounter the problem to reconcile the large values for T^ with 
the rather modest value for the resistivity. Using the leading order of the 1/N expansion 
it was shownil that the obtained values for Tc are comparable with those found in high- 
Tc superconductors. The question immediately arises whether the same theory (and also 
theories based on phonon- or spin fluctuation-mediated interactions) are in agreement with 
the observed low electronic scattering rates in these systems. For this aim we study the 
quantities 1/r and l/rtj. which describe the relaxation of quasi-particles and the momentum 
relaxation near the Fermi surface, the latter being intimately related to transport. Using 
our previously deflned effective particle-hole interaction G^'^ one obtainS 

l/r = — / deiV(e) «/meP^(k,k'-k,e) »o (nfe(e) + np(e)), (25) 

iV J —oo 



2 /-^^ , Jm0P'^(k,k'-k,e)(v(k) - v(k'))2 »o, , ^^ 

<< ... >>o denotes an average over the Fermi line with a constant eigenvector, i.e., the 
weight in the line integral around the Fermi line contains only the density of states. Im 
stands for the imaginary part, v(k) = 9e(k)/c}k is the velocity of the electron at the point 
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k. A^(e), nB^e) and npi^e) denote the density of states, the Bose and the Fermi distribution 
function, respectively. In the high temperature hmit 1/r and l/r^j. approach the hmits 2ttT\ 
and 27iT\tr, respectively, where A and Xtr are the corresponding dimensionless coupling 
strengths. In order to achieve the desired high values for with a spin-fluctuation or 
phonon mechanism A has to be roughly between two and three (2.6 in Ref.i, 3 in Ref.i^, 
2.28 using xp the values of Table II, second-last column, of Ref.i, 2.9 in Ref.il). On 
the other hand, infrared and resistivity data in many cuprates near optimal doping lead to 
a much lower experimental value for At,, of about 0.3 ± 0.1, see Table II in Ref.Ba. 

The dashed curves in Figs. 9 and 10 show 1/r and l/rtr, respectively, as a function of 
temperature for the kernel q'^Xp- The curves are quasilinear above 150 or 200 K correspond- 
ing to dimensionless coupling constants Ap and Xtr,p of about 2-3. Such a large value for Xp 
is needed in order to reproduce the large value for Tc. The corresponding transport value 
Xtr,p, however, is even larger by about 20 per cent. The reason for this is the following: If 
G is momentum independent A and Xtr are equal. Considering now kernels G which depend 
only on the difference of their momenta, i.e., k — k'. If G shows mainly forward scattering, 
i.e., if /mG(k — k') is large for small momentum and small at large momentum transfers 
Xtr is smaller than A. On the other hand if scattering occurs mainly with large momentum 
transfers A is smaller than Xtr- Spin-fluctuations are strongest in the region around Q in 
momentum space, i.e., at large momentum transfers. As a result Xtr^p > Ap in agreement 
with the slopes in Figs. 9 and 10. We expect that Xtr,p > Xp and A ~ 2 — 3 are generic 
values for spin fluctuation models for high-Tc superconductivity. Thus there is a severe 
discrepancy of about one order of magnitude between theoretical and experimental values 
for Xtr,p. 

The long-dashed curves in Figs. 9 and 10 describe relaxation rates due to phonons. 
Following our previous treatmentEy we used a Holstein model with a phonon frequency 
ujq = ASmeV and a dimensionless coupling strength Xph = 0.75 near optimal doping. As 
discussed in Refs.iil'0 vertex corrections must be included even in the leading order of the 
1/N expansion. The resulting vertex, taken at = 0, depends strongly on momentum in 
the region of the considered doping values. This implies that the constant, bare electro n- 
phonon coupling becomes effectively strongly momentum dependent, being large at small 
and small at large momentum transfers. As a result the phonon contribution to Atr^pjj. given 
essentially by the slope of the long-dashed line in Fig. 10, becomes much smallerE^ than 
Xph = 0.75, namely, Xtr,ph ~ 0.11. More generally, such a strong reduction in Xtr relative to 
A, should occur for all charge but not for spin interactions. This may explain why phonons 
and impurities play usually no role in transport data though their bare couplings to electrons 
may not be smallE3. 

The relation between Tc and relaxation rates is quite different in the case of the 1/N 
expansion. The instantaneous part of the kernel is large, for instance, its Fermi surface 
average in the d-wave channel is -0.48 at 5 = 0.17. This term alone would yield a of 
about the same magnitude as the whole 0(1) contribution for G. On the other hand, this 
instantaneous term does not contribute to the momentum relaxation , i.e., it does not enter 
transport. The solid lines in Figs. 9 and 10 represent the relaxation rates due to G^'' as given 
by Eqs.( pO| - p3D . The curves are again quasilinear at not too low temperatures. Their absolute 
values are rather small and correspond to dimensionless coupling constants of A ~ 0.25 and 
Xtr ~ 0.22. This means that charge and spin fluctuations do not scatter electrons very 
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efficiently near the Fermi surface. Finally, the long-dashed curves in Figs. 9 and 10 are the 
relaxation rates due to the spin fluctuation kernel 0^ of the 1 /N expansion. Comparing these 
curves to the dashed curves in Figs. 9 and 10 one recognizes that the leading spin fluctuation 
contribution, which is of 0(1/N) in B, yields only small scattering rates, compared to those 
in spin fluctuation models for high-Tc compounds. The total transport coupling constant is 
the sum of its electronic and phononic parts and about 0.5 in our model. This value is near 
the rather universal experimental value of 0.3 ± 0.10. 



VI. CONCLUSIONS 

Using a 1/A^ expansion for the t — J model in terms of X operators the 0(1/N) contri- 
butions to the anomalous self-energy just below and the self-energy in the normal state 
have been considered and numerically evaluated for typical values of the model parameters. 
Our treatment is restricted to the optimally doped and the overdoped region defined by the 
absence of additional instabilities, e.g., towards a flux state. Unlike to the case of phonons 
or spin fluctuations the resulting effective interactions in the particle-particle and particle- 
hole channels turn out to be different from each other. We also show how scattering rates 
relevant to transport can be obtained form the expression for the self-energy. In order to 
examine the convergence of the 1/A^ expansion and to take into account antiferromagnetic 
effects we also have included graphs of 0{1/N'^) which describe spin-fluctuation mediated 
contributions to the effective interaction between electrons. Though these terms are only a 
subset of all l/N"^ contributions it seems not unreasonable to assume that they represent 
the leading corrections to the 0(1/N) terms if, as is often believed, spin fluctuations play 
an important role in the t — J model in the considered doping region. 

Our main conclusion is that the leading 0(1/N^) contribution describing spin fluctuations 
gives only a small correction to the effective interaction and obtained already in the leading 
order 0(1/N). As discussed in Ref.i the effective interaction in 0(1/N) contains no RPA-like 
spin fluctuations. Instead, this effective interaction and thus also d-wave pairing is caused 
by the instantaneous contribution of the Heisenberg term and by spin flip processes induced 
by the X operator algebra and thus by the constraint. The reason for the smallness of the 
0(1/N^) contribution is that the momentum averaged magnetic susceptibility at large N's, 
Xoo, is, as a function of frequency, rather structureless and extends over the whole bandwidth. 
Though the effective coupling constant between electrons and spin excitations is set by t 
( and not J) the resulting contribution to the imaginary part of the effective interaction 
(the analogue of a^F of the phonon case) is small compared to the spin flip (and charge) 
contribution obtained in 0(1/N). In contrast to that the spin susceptibilities used in spin 
fluctuations theories such as xp of Refs.i^ is large and strongly peaked at low frequencies. 
In other words, most of the spectral weight is in this case accumulated at small frequencies 
of the order of a fraction of J whereas it is rather constant and spread over an energy scale 
of t in the case of Xoo- Available magnetic susceptibility spectra from inelastic neutron 
scattering in Y Ba2Cu^0'!-x refer mostly to the strongly underdoped region and resemble 
Xp more than Xoo- In optimally doped samples of Y Ba2Cu^0-j-x spin scattering is too weak 
to be detected presently. The obtained experimental upper bound for x", however, suggests 
in view of the sum rule that x" is rather structureless and spreads over an energy scale set 
by t similar as x'Lo does. There is another fundamental difference between calculations based 
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on 1/N expansions and spin fluctuation models: The effective interaction of 0(1/N) in the 
d-wavc channel contains strongly pair-breaking charge excitations at high energies of the 
order t which cancel the contribution of pair-forming low-energy spin flip excitations to a 
large extent. These pair-breaking terms are entirely neglected in the spin fluctuation models 
but play an important role for quantitative calculations of Tc in the 1 /N expansion. 

For the same reasons as above wc find that the leading spin fiuctuation terms of 0(1/N^) 
modify only little the scattering rates 1/r and l/r^r obtained already in 0(1/N). If spin fiuc- 
tuation terms would dominate and cause the large values of Tc the corresponding scattering 
rate l/rtr would exceed the experimental ones by far. Responsible for this is also the fact 
that the spin susceptibility is large for large momentum transfers which makes l/rtr even 
larger than 1/r. In the 1/N expansion large values for Tc are compatible with small values 
for l/rtr'- Spin fiuctuation terms are small and Tc is to a large extent determined by the 
instantaneous part of the interaction which does not contribute to relaxation rates. The case 
of phonons or, more generally, of charge-charge interactions is somewhat different from spin 
interactions. Here the charge vertex due to correlation effects has to be included already 
in leading order of the 1/N expansion. As a result collective effects dominate and charge 
interactions are confined to small momenta causing l/rtr << 1/r near optimal doping. This 
may explain why phonons and impurities play only a minor role in the transport data of 
optimally doped high-Tc compounds. 
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FIGURE CAPTIONS 



Fig.l: Momentum integrated magnetic susceptibility 'x"{uj) in units of n^/eV. Dot-dashed 
line: t — J model for — * oo; dashed line: ^'p{u) of Refi; circles: exp. points for 
YBa2Cu^0Q,Q at T = SOi^B. 

Fig. 2: Magnetic susceptibility x"(Q,co') for Q = (vr, vr) in units of ii\/eV. Dot-dashed line: 
t — J model for — > cxD multiplied by 10; dashed line: Xp(Q)'^) of Ref.i; circles: exp. points 
for YBa2Cu^0e,^ at T = 200idll. 

Fig. 3: Absolute value of the real part of the susceptibility at = versus distance. Circles 
and diamonds correspond to Xoo and xpi respectively, at neighboring sites of the square 
lattice, the dashed and dot-dashed lines to Kq{x/^) with ^ = 2.3 and ^ = 1.0, respectively. 

Fig. 4: Dependence of the static kernel 0(k, k',0) on the second momentum k' along the 
Fermi line for a fixed first momentum k = (2.465,0.309); solid line: 0(1) contrib.; dot- 
dashed line: 0(1/N) contrib.; dashed line: g^XP divided by 8. 

Fig. 5: Lowest eigenvalue of static kernels in (d-wave) symmetry; solid, dashed, dot- 
dashed, and dotted lines correspond to 0(1) of 9, g^Xp^ the numerator, and denominator 
of Eq.(18), respectively. The dot-dashed curve corresponds at the same time to 0^. 

Fig. 6: Negative imaginary part of the d-wave projected kernel Qz{uj + irf) as a function of 
the frequency; sohd, dot-dashed, and dashed lines correspond to the contributions of 0(1), 
0(1/N), and q^Xp to 63, respectively. 

Fig. 7: Transition temperature in Kelvin as a function of doping; circles, crosses, and 
diamonds correspond to the contributions of 0(1), 0(1/N), and 0(1)+0(1/N), respectively. 

Fig. 8: Negative imaginary part of the 0(1/N) electronic self-energy as a function of energy 
with the momentum averaged over the Fermi surface. Dot-dashed, dashed, and solid lines 
describe the charge, spin and total contributions, respectively. 

Fig. 9: Quasi-particle relaxation rate 1/r as a function of temperature. Dashed, long-dashed, 
solid, and dot-dashed lines correspond to xpi an electron-phonon model, and the 0(1/N) 
and 0(1/N^) contributions to the self-energy of the t — J model, respectively. 

Fig. 10: Transport relaxation rate l/rtr as a function of temperature. Dashed, long-dashed, 
solid, and dot-dashed lines have been calculated from xpi an electron-phonon model, and 
the 0(1/N) and 0(1/N^) contributions to the self-energy of the t — J model, respectively. 
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